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<D . Abstract 

Fruitful ideas on how to quantize gravity are few and far between. In this 
paper, we give a complete description of a recently introduced non-perturbative 
gravitational path integral whose continuum limit has already been investigated 
extensively in d < 4, with promising results. It is based on a simplicial regulari- 
zation of Lorentzian space-times and, most importantly, possesses a well-defined, 
non-perturbative Wick rotation. We present a detailed analysis of the geometric 
and mathematical properties of the discretized model in d = 3, 4. This includes 
a derivation of Lorentzian simplicial manifold constraints, the gravitational actions 
and their Wick rotation. We define a transfer matrix for the system and show that 
it leads to a well-defined self-adjoint Hamiltonian. In view of numerical simulations, 
we also suggest sets of Lorentzian Monte Carlo moves. We demonstrate that certain 
pathological phases found previously in Euclidean models of dynamical triangula- 
tions cannot be realized in the Lorentzian case. 
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1 Introduction 

Despite an ever-increasing arsenal of sophisticated mathematical machinery, the non- 
perturbative quantization of gravity remains an elusive goal to theorists. Although 
seemingly negligible in most physical situations, quantum-gravitational phenomena 
may well provide a key to a more profound understanding of nature. Unfortunately, 
the entanglement of technical problems with more fundamental issues concerning 
the structure of a theory of quantum gravity often makes it difficult to pinpoint why 
any particular quantization program has not been successful. 

The failure of perturbative methods to define a fundamental theory that includes 
gravity has led to alternative, non-perturbative approaches which seek to describe 
the quantum dynamics of gravity on the mother of all spaces, the "space of geome- 
tries". By a "geometry" we mean a space-time with Lorentzian metric properties. 
Classically, such space-times usually come as smooth manifolds M equipped with a 
metric tensor field g^ v {x). Since any two such metrics are physically equivalent if 
they can be mapped onto each other by a diffeomorphism of M, the physical degrees 
of freedom are precisely the "geometries" , namely, the equivalence classes of metrics 
with respect to the action of Diff(M). 

In what follows, we will describe a path-integral approach to quantum gravity 
that works much more directly with the "geometries" themselves, without using 
coordinates^. It is related in spirit to other constructions based on simplicial "Regge" 
geometries ||, in particular, the method of dynamical triangulations j|, |5| (see |J 
for a comparison and critical appraisal of various discrete approaches to 4d quantum 
gravity). Its main aim is to construct an interacting quantum theory of gravity in 
four dimensions as the continuum limit of a well-defined statistical mechanics model 
of regularized gravitational field configurations. The formal continuum path integral 
for gravity is represented by a discrete sum over inequivalent triangulations T, 



where m(T) is a measure on the space of discrete geometries. Unlike previous 
approaches, we use a space of piecewise linear Lorentzian space-times as our starting 
point. That this procedure is in general inequivalent to path integrals over bona 
fide Euclidean geometries has already been demonstrated in two dimensions [|7|, ||. 

The reason why previous discrete path-integral approaches have been formu- 
lated for Euclidean gravity^] is often a technical one, rather than the conviction that 
Euclidean space-times are more fundamental physical quantities than Lorentzian 

2 There is no residual gauge invariance since the state sum is taken over inequivalent discretized 
geometries. In this sense, the formalism is manifestly diffeomorphism-invariant. In Euclidean 
2d quantum gravity, there is ample evidence that this procedure is equivalent to a gauge-fixed 
continuum formulation, see ||, |^] for reviews. 

3 Euclidean gravity is the theory with (minus the) action (^|) below, where the Lorentz metric 
g^y has been substituted by a Riemannian metric <?^" cl with positive-definite signature. 
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ones (after all, a phenomenon like black holes is intimately related to the causal, 
Lorentzian structure). The corresponding real weight factors e~ scuc are then used 
in Monte Carlo simulations, and also usually necessary (though not sufficient) for a 
convergence of the non-perturbative state sums. 

The problem with nonperturbative Euclidean quantum gravity is two-fold. First, 
to our knowledge and aside from our own proposal, no Wick rotation on the space of 
all geometries has ever been constructed in a path- integral context, leaving unclear 
its relation with the physical theory of Lorentzian gravity. Although it is logically 
possible that a natural Wick rotation may suggest itself once the Euclidean quantum 
theory has been solved, there is so far no convincing evidence (in dimension four) 
that a non-trivial continuum theory exists. This is the second problem of Euclidean 
quantum gravity. These conclusions seem to hold for a wide variety of approaches, 
including those based on gauge-theoretic instead of geometric variables |J. 

Our alternative proposal consists in taking the Lorentzian structure seriously 
from the outset, using a well-defined Wick rotation on the full, regularized path 
integral ([!]), performing all calculations, simulations and continuum limits in the 
"Euclidean sector" , and finally rotating back the results. As mentioned above, this 
program has been carried out successfully and completely explicitly in d = 2 J/J. 
Another lesson from two-dimensional quantum gravity is that the Lorentzian and 
Euclidean theories will in general be inequivalent, even if one can Wick-rotate the 
final propagator of the Euclidean theory, in this case, of the Euclidean Liouville 
theory (see || for an overview of 2d Lorentzian results). Although we do not wish 
to exclude the possibility that interesting theories of Euclidean quantum gravity in 
higher dimensions exist, there seems to be a shortage of ideas of how to modify 
current models to lead to an interesting continuum behaviour |H| . 



The apparent inequivalence of our Lorentzian proposal makes it a promising and 
genuinely new candidate for a theory of quantum gravity. More concretely, the 
geometric constraints which are a natural consequence of the causal structure act 
as a "regulator" for the path integral, suppressing the types of highly degenerate 
geometries that seem to be related to the failure to date of the Euclidean approaches, 
be they gauge-theoretic |yj , of Regge type |12) , or dynamically triangulated [K| . 



to produce any convincing evidence of interesting continuum physicsQ. This issue is 
discussed further in Sec. |8| below. 

In our construction, we mostly use conventional tools from quantum field theory 
and the theory of critical phenomena applied to a diffeomorphism-invariant theory. 
However, it does not take much to realize that general coordinate invariance cannot 
be accomodated by merely making minor changes in the standard QFT formalism. 
Almost all methods of regularizing and renormalizing make use of a fixed metric 
background structure (usually that of flat Minkowski space), albeit often in an im- 
plicit way. By contrast, in a non-perturbative formulation of quantum gravity the 



4 In Euclidean quantum Regge calculus, an unusual prescription has recently been s ugg ested 
to work around the problem of degenerate geometries at weak gravitational coupling [[L4|: the 
continuum limit is to be reached by an analytic continuation from the strong-coupling phase. 
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metric degrees of freedom are dynamical and initially all on the same footing. A 
"ground state of geometry" should only emerge as a solution to the quantum equa- 
tions of motion, much in the same way as a classical physical space-time is the result 
of solving the classical Einstein equations. 

As indicated in ([]]), our path integral will be given as a sum over discretized 
Lorentzian geometries (for given boundary data), each weighed with e lS , where S 
is the Regge version of the classical Lorentzian gravitational action in d space-time 
dimensions (d = 3,4), 

S[ g ^] = ^ f d d x yj- det g(R — 2A) + k [ d d ~ x x Vdet hK, (2) 

2 J M JdM 

with k" 1 = 87rGNewton and the cosmological constant A. We have included a stan- 
dard boundary term depending on the induced metric hij on the boundary and its 
extrinsic curvature Ky. In order to make the state sum well-defined, we use a non- 
perturbative Wick rotation on the discretized geometries to convert the oscillating 
amplitudes into real Boltzmann factors e~ 5 ° uc . It can then be shown that the en- 
tire non-perturbative state sum converges for suitable choices of the bare coupling 
constants, even in the infinite- volume limit. 

It should be emphasized that our construction is intrinsically quantum. The 
formalism of dynamical triangulations is not well suited for approximating smooth 
classical space-times, because it lacks the notion of an infinitesimal smooth variation 
of its field variables^. This would be required if we wanted to evolve the classical 
Einstein equations numerically, say. It by no means disqualifies the method from 
use in the quantum theory. There, one is faced with the problem of having to ap- 
proximate "ergodically" the space of all geometries, whereas in classical applications 
one is usually interested in approximating individual, smooth geometries (often so- 
lutions to the equations of motion). One can compare this situation to the Feynman 
path integral for a non-relativistic particle, which can be obtained as the continuum 
limit of a sum over piecewise linear trajectories. Although the final state sum is 
dominated by nowhere differentiable (and therefore highly non-classical) paths, this 
is no obstacle to retrieving (semi-) classical information about the system from the 
expectation values of suitable observables. 

Most of the results presented below were already announced in [fT5| , where also 



some of the 3d formulas were given explicitly. Our account here is comprehensive 
in that we give all relevant calculations in both three and four dimensions. 

We start by a brief description of how to obtain Lorentzian space-times by gluing 
together flat simplicial Minkowskian building blocks. We compute their volumes and 
dihedral angles as functions of a, where —a denotes the squared (geodesic) length 
of one- dimensional time-like edges, and the space-like edge lengths are fixed to 1. 
In the next section, we derive some topological identities among the bulk numbers 



5 This is different from the situation in (quantum) Regge calculus, where continuous changes in 
the edge length variables are allowed. 
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Ni, i = 1, . . . , d of simplices of dimension i. These generalize the Euler and Dehn- 
Sommerville identities to the Lorentzian situation where we distinguish between the 
space- and time-like character of the building blocks. These identities are used in 
the next section, where the gravitational actions are computed as functions of a, 
and where we demonstrate how the corresponding Euclidean actions are obtained 
by a suitable analytic continuation of a to negative real values. 

With these ingredients in hand, we define the transfer matrix of the regularized 
model in Sec. [|, together with the Hilbert space it acts on. Some important prop- 
erties of the transfer matrix are proved in Sec. || which ensure the existence of a 
well-defined Hamiltonian. Paving the way for numerical simulations of the model, we 
discuss possible sets of geometric Monte Carlo moves in Sec. |7|. Using these moves, 
one can understand certain "extreme" regions of the kinematical phase space. We 
show in Sec. |] that the regions of highly degenerate geometries of Euclidean quan- 
tum gravity cannot be reached by Lorentzian dynamical triangulations. We end 
with a summary of current achievements and an outlook on the potential applica- 
tions of discrete Lorentzian gravity. The two appendices contain a brief reminder of 
some properties of Lorentzian angles and a proof of link-reflection positivity of 2d 
Lorentzian gravity. 



2 Discrete Lorentzian space-times 

Our first task will be to define which discrete Lorentzian geometries T contribute to 
the sum in the path integral (|l|). They are a straightforward generalization of the 



space-time histories used in two dimensions, and have already been described [|15 
and used |16|] in three dimensions. 



Briefly, they can be characterized as "globally hyperbolic" d- dimensional sim- 
plicial manifolds with a sliced structure, where (d — l)-dimensional "spatial hyper- 
surfaces" (i.e. Euclidean (d— l)-dimensional equilaterally triangulated manifolds) 
of fixed topology are connected by suitable sets of d- dimensional simplices. As a 
concession to causality, we do not allow the spatial slices to change topology. There 
is a preferred notion of a discrete "time" , namely, the parameter labelling successive 
spatial slices. Note that this has nothing to do with a gauge choice, since we were 
not using coordinates in the first place. This "proper time" is simply part of the 
invariant geometric data common to each of the Lorentzian geometries. 

We choose particular sets of simplicial building blocks in three and four dimen- 
sions. Unlike in the Euclidean case, we cannot make the simplest possible choice of 
making all <i-simplices equilateral. Instead, we fix all spatial (squared) link lengths 
to 1, ^pace = 1; an d all time-like links to have a squared length /^ ime = —a, a > 0. 
Keeping a variable allows for a relative scaling of space- and time-like lengths and 
is convenient when discussing the Wick rotation later. As usual, the simplices are 
taken to be pieces of flat Minkowski space, and a simplicial manifold acquires non- 
trivial curvature through the way the individual building blocks are glued together. 
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For simplicity, we have set the "lattice spacing" a, which defines a diffeomorphism- 
invariant cutoff for our model, equal to 1. (In general, we should have written 
/g pace = a 2 and l^ ime = —aa 2 . The correct dependence on a in the formulae below can 
be inferred from dimensional considerations.) 

As usual in the study of critical phenomena, we expect the final continuum the- 
ory (if it exists) to be largely independent of the details of the chosen discretization^. 
Our choice of building blocks is simple and allows for a straightforward Wick rota- 
tion. However, other types of fundamental building blocks may sometimes be more 



convenient. For example, a combination of pyramids and tetrahedra is used in [[IS 
in a three-dimensional context. In general, one should be careful not to restrict the 
local curvature degrees of freedom too much. With our choice of building blocks, 
the (discretized) local curvatures around given bones (simplices of dimension d—2) 
can always take both positive and negative values. 

We will now go on to compute the volumes and dihedral angles of the ci-dimen- 
sional Minkowskian simplices, because they are needed in the gravitational Regge 
action in dimension d. All (Lorentzian) volumes in any dimension we will be using 
are by definition real and positive. Formulas for Euclidean volumes and dihedral 
angles can be derived from elementary geometric arguments and may be found in 
many places in the literature p0|. They may be continued to Lorentzian geometries 



by taking suitable care of factors of i and —1. We will follow Sorkin's treatment 



and conventions for the Lorentzian case (Some basic facts about Lorentzian 



angles are summarized in Appendix 1). The dihedral angles B are chosen such that 
< Re0 < 7r, so giving sin0 and cosB fixes them uniquely. The angles are in 
general complex, but everything can be arranged so that the action comes out real 
in the end, as we shall see. 

A further ingredient needed in the calculations are the volumes of lower-dimen- 
sional simplices. By convention, we have Vol(point)= 1. For the one-dimensional 
links, 

Vol(space-like link) = 1 and Vol(time-like link) = y/a. (3) 

Also for the triangles, we must distinguish between space- and time-like. The former 
lie entirely in planes t = const whereas the latter extrapolate between two such slices. 
Their respective volumes are 

Vol (space-like triangle) = ^ and Vol(time-like triangle) = jy/Aa + 1. (4) 
2.1 Geometry of three-simplices 

The simplices of top-dimension are tetrahedra, of which (up to reflection symmetry) 
there are two types, (3,1) and (2,2). Starting with the former (Fig. la), one computes 



Vol(3, 1) = — VSgTTT. (5) 



6 In dimension two, this has already been shown for both the Euclidean |17 and the Lorentzian 
flStl theories. 
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^ (a) (b) 

Figure 1: A (3,1)- and a (2,2)-tetrahedron in three dimensions. 
The dihedral angle around any of the space-like bones is given by 



cos 9( "» = Sm e(s '" = VJM&TT <6) 



and around a time-like bone by 



n 2a + 1 • n gyW3Q + 1 

cos (3.D = I^TT' Sm e ^> = 4a + l • (7) 



For the (2,2)-tetrahedra (Fig. lb) one finds 

Vol(2,2) = ^i=v / 2^+T. (8) 
The dihedral angle around a space-like bone is fixed by 

n 4a + 3 • n • ^V^oT+T 

{ ' > 4a + 1' ( , ) 4a + 1 ' K ' 

and that around a time-like bone by 



cos e (2>2) = sm e (2 , 2) = 4a + 1 • (io) 

2.2 Geometry of four-simplices 

In d — 4 there are up to reflection symmetry two types of four-simplices, (4,1) and 
(3,2). In addition to the volumes computed above, we now also need the volumes 
of equilateral tetrahedra that lie entirely in slices t = const, given by Vol (space- 
like tetrahedron)= For the volume of the simplices of type (4,1), Fig. 2a, one 
computes 

Vol(4,l) = ^ v / 8oT3. (11) 

A (4,l)-simplex contributes to the curvature with two types of dihedral angles 
(of which there are a total of 10 per simplex), depending on whether the parallel 
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Figure 2: A (4,1)- and a (3,2)-tetrahedron in four dimensions. 



transport is taken around a space- or a time-like triangle. There are six of the 
former, each contributing a dihedral angle determined by 



cos G( 4i i) = -—=—===, sin 6(4,1) = \ j ofn _ , ^ , (12) 



'3(8a + 3) 

2 v / 2\/3^TT y 8(3a + 1 

and similarly for the angle around time-like triangles 



n 2a + 1 . V4a + lV8a + 3 

COSC>C4il = — r, Sill 9(4 n = ; . (13) 

2(3a + l)' ( '> 2(3a + l) 

As in the three-dimensional case (f|) above, the trigonometric functions ([l|) are 
complex, with a constant real part equal to tt/2, whereas (|13"D for the Euclidean 
angles are real. 

The other type of four-simplex is the (3,2)-simplex, see Fig. 2b. Its volume is 
given by 

Vol(3,2) = -^Vl2a + 7. (14) 

As for the dihedral angles, the situation is slightly more involved. Of the ten angles 
contributing per four-simplex there is exactly one around a space-like triangle, with 

6a + 5 ^ ^3^l2a + 7 

o~7q — rTT' sm 9(3,2) = -% -rrr^ — —rr 
2{3a + 1) 2(3a + 1) 

which are purely imaginary. There are two varieties of angles around time-like 
triangles. If the angle is formed by two three-dimensional "faces" that are both 
tetrahedra of type (2,2), one has 

Aa + 3 . V(4a + l)(12a + 7) 
C ° S @(3 ' 2) = 4(2^+1) ' Sm @(3 ' 2) = 4(to + l) • (16) 

By contrast, if the two tetrahedra involved are a pair of a (3,1) and a (2,2), the 
angle between them is defined by 



COS 0(3,2) = — r — , Sin 0(3,2) = ~% 0/0 „ , 1N , (15) 



-1 n J 4a + 1 12a + 7 

cos ( 3 2 ) = — , sin 0o 2) = v , , , 17 

2 v / 2v / 2an v / 3^Tl { ' } 2v / 2v / 2an v / 3^TT 
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3 Topological identities for Lorentzian triangula- 
tions 



In this section we derive some important linear relations among the "bulk" variables 
Ni, i = 0, . . . , d which count the numbers of 2-dimensional simplices in a given d- 
dimensional Lorentzian triangulation. Such identities are familiar from Euclidean 
dynamically triangulated manifolds (see, for example, [22|, ^3|). The best-known of 
them is the Euler identity 

x = Nq - Nl + N 2 - N 3 + . . . , (is) 

for the Euler characteristic x of a simplicial manifold with or without boundary. For 
our purposes, we will need refined versions where the simplices are distinguished by 
their Lorentzian properties. The origin of these relations lies in the simplicial mani- 
fold structure. They can be derived in a systematic way by establishing relations 
among simplicial building blocks in local neighbourhoods and by summing them 
over the entire triangulation. Our notation for the numbers Ni is 

Nq = number of vertices 
N± L = number of time-like links 
N± L = number of space-like links 
Nj L = number of time-like triangles 
iVf L = number of space-like triangles 

Nj Ll = iVg 3 ' 1 '' = number of time-like (3,1)- and (l,3)-tetrahedra 
Nj h2 = iVg 2 ' 2 '' = number of time-like (2,2)-tetrahedra 
iVf L = number of space-like tetrahedra 

Nj Ll = iVj = number of time-like (4,1)- and (l,4)-simplices 

Nj L ' 2 = A^ 3 ' 2 " 1 = number of time-like (3,2)- and (2,3)-simplices. (19) 



3.1 Identities in 2+1 dimensions 

We will be considering compact spatial slices ( 2 )£, and either open or periodic bound- 
ary conditions in time-direction. The relevant space-time topologies are therefore 
J x ' 2 )£ (with an initial and a final spatial surface) and S 1 x^E. Since the latter 
results in a closed three-manifold, its Euler characteristic vanishes. From this we 
derive immediately that 

x(Ix {2) Z)=x( {2) Z). (20) 

(Recall also that for closed two- manifolds with g handles, we have x = 2 — 2g, for 
example, x(S 2 ) = 2 for the two- sphere.) 

Let us for simplicity consider the case of periodic boundary conditions. A three- 
dimensional closed triangulation is characterized by the seven numbers Nq, Nf L , 
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N? L , Nl h , Nj L , N^' 1] and ivf ,2) . Two relat ions among them are directly inherited 
from the Euclidean case, namely, 

iVo - iVf L - N? L + iV 2 SL + Nj L - N^' 1] - ivf ,2) = 0, (21) 
iV 2 SL + Nj L = 2(iV 3 (3 ' 1) + iV 3 (2 ' 2) ). (22) 

Next, since each space-like triangle is shared by two (3,l)-tetrahedra, we have 

jvf' 1 ) = ijvf\ (23) 
3 

Lastly, from identities satisfied by the two-dimensional spatial slices, one derives 

iVf = ~N$ h , (24) 

N = X ^Z)t+±Nl L , (25) 

where we have introduced the notation t for the number of time-slices in the trian- 
gulation. 

We therefore have five linearly independent conditions on the seven variables iVj, 
leaving us with two "bulk" degrees of freedom, a situation identical to the case of 
Euclidean dynamical triangulations. (The variable t does not have the same status 
as the Ni, since it scales (canonically) only like a length, and not like a volume.) 



3.2 Identities in 3+1 dimensions 

Here we are interested in four-manifolds which are of the form of a product of a 
compact three- manifold ^E with either an open interval or a circle, that is I x ®E 
or S 1 x ( 3 )E. (Note that because of x( (3) S) = 0, we have x(I x (3) S) = x(^ 1 x (3) S). 
An example is xi^ 1 x T 3 ) = x{T 4 ) = 0.) In four dimensions, we need the entire 
set ( ]T9| ) of ten bulk variables iVj. Let us again discuss the linear constraints among 
them for the case of periodic boundary conditions in time. 

There are three constraints which are inherited from the Dehn-Sommerville con- 
ditions for general four-dimensional triangulations fl22| , ^3| , 

N - Nf L - N? L + Nl L + Nj L - iV 3 SL - iV 3 TLl - Nj L2 + iV 4 TLl + A^ 2 = x , 
2(A 1 SL + A 1 TL ) - 3(A 2 SL + A 2 TL ) + A(Nl h + Nj Ll +Nj L2 ) - 5(Nj Ll +Nj u ) = 0, 
5(A 4 TLl + Nj u ) = 2(A 3 SL + Nj u + Nj u ). (26) 

The remaining constraints are special to the sliced, Lorentzian space-times we are 
using. There are two which arise from conditions on the space-like geometries alone 
(cf. (0), (^)), 

N - Nf L + iV 2 SL - iVf L = 0, 

JV 2 SL = 2iV 3 SL . (27) 
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Furthermore, since each space-like tetrahedron is shared by a pair of a (4,1)- and a 
(l,4)-simplex, 

2jVf L = iV 4 (4,1) , (28) 

and since each time-like tetrahedron of type (2,2) is shared by a pair of (3,2)- 
simplices, we have 

2Nj L2 = 3ivf' 2) . (29) 

In total, these are seven constraints for ten variables. This is to be contrasted 
with the case of Euclidean triangulations, where there are only two bulk variables. 
With the help of the results of Sec. |8| below, one can convince oneself that there are 
no further linear constraints in the Lorentzian case. Taking iV 2 TL , jvf ' 1} and ivf ' 2) 
as the three remaining variables, it is straightforward to verify that a relation of the 
form 

aiV 2 TL + (3Njf' 1] + 1 Nf' 2) +5 = (30) 
is not compatible with the Monte Carlo moves, (|S2j)-(|8~5|), unless a = (3 = ^ = 5 = 0. 



4 Actions and the Wick rotation 

We are now ready to construct the gravitational actions of Lorentzian dynamical 
triangulations explicitly. We will investigate their behaviour in the complex a-plane 
and show that - subject to a dimension-dependent lower bound on a - the map 

1 7 

a i— > —a, a > — (d = 3), a > — (d = 4), (31) 

where a is continued through the lower half of the complex plane, defines a non- 
perturbative Wick rotation from Lorentzian to Euclidean discrete geometries. Under 
this map, the weights in the partition function transform according to 



e lS ^e~ scucl . (32) 

For the special case a = — 1, the usual expressions for the actions used in equilateral 
Euclidean dynamical triangulations are reproduced. 

4.1 Gravitational action in three dimensions 

The Regge analogue of the continuum action (g) in d = 3 is given by (c.f. |2"T||) 



= k Vo1(0t(2tt- ) + A: Vo1(0(2tt- ^ 6) 

space — like tetrahcdra time — like tetrahcdra 

I at / / at I 

-A Vol(3,l)-A Vol (2,2), (33) 



(3,1)&(1,3)- (2,2)- 
tetrahedra tetrahedra 
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where we have now set A = kA. Performing the sums, and taking into account how 
many tetrahedra meet at the individual links, one can re-express the action as a 
function of the bulk variables iV\ and N 3 , namely, 

= k ( ^iVf L - 2 -Nf» arcsin 2 ^^+ l _ ^(M) arccog . 



i i 4a + 1 z y/Zy/Aa + 1 

( ~ iN? * 2) arccos - 3ivf ' 1} arccos ^±1) 

V 4a + 1 d 4a + 1 / 

-A (ivf ' 2) ±V^+2 + ^v/SaTl) • (34) 

Our choice for the inverse trigonometric functions with imaginary argument from 
®' ®> av °ids branch-cut ambiguities for real, positive a. Despite its appearance, 
the action (|34]) is rea/ in the relevant range a > 0, as can be seen by applying 
elementary trigonometric identities and the relation The final result for the 

Lorentzian action can be written as a function of three bulk variables (c.f. Sec. |3|), 
for example, N]- L , and N^ 2 ' 2 \ as 

_l_ jy^ 3 ' 1 ) arcsinh — — 3k\/a arccos - — — V3a + 1 

3 V V^V^a + 1 4a + 1 12 

+ ivi 2 ' 2 "* (2k arcsinh — ~ ^ a ^ — Ak\fa. arccos — V4a + 2 

3 V 4a + 1 v 4a + 1 12 

(35) 

Since we are interested in continuing this action to the Euclidean sector where 
a is real and negative, we need to understand the properties of fl3*5D in the complex 
a-plane. For general a (and assuming k and A real and arbitrary), S^(a) takes 
complex values. There are two special cases: (i) for a G H, a > 0, the entire action 
is real - this is simply the Lorentzian case considered above; (ii) for a E M, a < — |, 
the expression (|35D is purely imaginary and we will see that it coincides with iS end . 
The reason for the non-trivial upper bound of — 5 are the triangle inequalities for 
Euclidean signature. It is easy to verify that as the positive Euclidean squared 
length —a becomes shorter, the (2,2)-tetrahedral building blocks degenerate first, 
at —a = |, followed by the (3,l)-tetrahedra at —a = |. Therefore, if for fixed a we 
adopt a 1 — > —a as prescription for a Wick rotation, we should confine ourselves to 
values |a| > |. 

Because of the appearance of square roots in the action Q3"5|), there is a two-fold 
ambiguity in what we mean by S^^—a). It turns out that the correct Euclidean 
expression is obtained by continuing a through the lower-half complex plane such 
that the branch cuts along the negative real axis are approached from below. In 
other words, if the argument of an expression of the form y/a + c becomes negative 
as a result of the Wick rotation, we interpret it as — i\J — (a + c). 
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Let us now consider the special case a = 1, together with its Wick rotation. The 
Lorentzian action becomes 

S^\ a = 1) = 2nkNj h + (-3.548k - 0.167A) + JVf ' 2) (-5.355A: - 0.204A) 
= 2irk X T + N^' 1] (-0A07k - 0.167A) + ivf ' 2) (0.929A; - 0.204A), 

(36) 

where in the last step we have used the Lorentzian manifold identities derived in 
Sec. ||. Its Wick-rotated Euclidean counterpart is given by 

S®(a =-1) = -2mkN? h + N^ 1] (2.673ifc + 0. 1 18zA) + ivf' 2) (7.386ifc + 0.1 18zA) 
= -2mk(N? h + Nf L ) + (A^3 (3,1) + iV 3 (2 ' 2) )(7.386iA; + O.II82A) 

= -2mkNi + N 3 (6ik arccos - + — ^A) = iS^ T , (37) 

3 6y2 

where in the last step we have identified the well-known expression for the action of 
Euclidean dynamically triangulated gravity in three dimensions. 

4.2 Gravitational action in four dimensions 

The form of the discrete action in four dimensions is completely analogous to (|33|) , 
that is, 

= k Vo1(A)-(2tt- J2 + k Vo1(A)(2tt- ) 

space — like 4 — simpliccs time — like 4 — simpliccs 

A at A A at A 

-A J2 Vol(4,l)-A J2 Vol (3,2). (38) 

(4,1)&(1,4)- (3,2)&(2,3)- 
tetrahedra tetrahedra 

Expressed in terms of the bulk variables N2 and N±, the action reads 

S« = k(^N*> - arcsin ~* ^ 12a + 7 

V i 4 2 U 4 2(3a+ 1) 



-^N^ arccos ) + -y^Tl ( 2nNj L 

AT^ 2 )(a ~ l 1 4a + 3 

—A/4 I o arccos - — — + 3 arccos ■ 



2 v / 2v / 2« + V30 + 1 4(2a + 1) 



r( 4 .i) Q _ c 2o+l \ / ( 4,i) ^80 + 3 (3) 2) Vl2a + 7 



We have again taken care in choosing the inverse functions of the Lorentzian angles 
in ([12]) and ([15]) that make the expression (|39|) unambiguous. Using the manifold 
identities for four- dimensional simplicial Lorentzian triangulations derived in Sec. H, 



13 



the action can be rewritten as a function of the three bulk variables Nj L , AQ ' and 
, in a way that makes its real nature explicit, 



2 



(4,1) 
4 



(-V3k arcsinh * - ^v^TTarccos ** + \ - A^±5) 

V 2v / 2v / 3^TI 2 V 2(3a + 1) 96 J 



V3k . . \/3Vl2a + 7 3A; 



-1 4a + 3 \ Jl2a + 7\ , . 

2 arccos — - + arccos — — A . 40 

2 v / 2 v / 2^Tl\/3a T TI 4(2a + l)7 96 J V ; 

It is straightforward to verify that this action is real for real a > — |, and purely 
imaginary for a G H, a < — -h. Note that this implies that we could in the 
Lorentzian case choose to work with building blocks possessing light-like (null) edges 
(a = 0) instead of time-like edges, or even work entirely with building blocks whose 
edges are all space-like. However, if we want to stick to our simple prescription 



for the Wick rotation, we should choose only values \a\ > Only in this case will 
we avoid degeneracies in the Euclidean regime in the form of violations of triangle 
inequalities. (Coming from large positive length-squared —a, it is the (3,2)-simplices 
that degenerate first, at —a = ^, followed by the (4,l)-simplices at —a = |.) 

As in three dimensions, in order to obtain the Euclidean version of the action, 
we continue the square roots in fl40|) in the lower half of the complex a-plane. Let 
us again consider the two special cases a = ±1. Evaluating flip]) numerically for the 
Lorentzian case, we obtain 

S( 4 ) ( a = 1) = 3.512kNj L + N^ A) (-4.284A; - 0.035A) + A"f ' 2) (-6.837A; - 0.045A) 
= 7.025^ + A"i 4,1) (-0.772A; - 0.035A) + n[ 3 ' 2) (0.188A; - 0.045A), 

(41) 

where we have used some topological identities from Sec. ||]. Its Wick-rotated version 
is given by 

S {A) (a = -l) = -2.721ikNj L + N^ 1] (2.987ik + 0.023i\) + 
A^f* 2) (5.708iA; + 0.023iA) 
= -2.721ifc(A" 2 TL + Nl L ) + (Af' 1} + A"f ' 2) )(5.708^ + 0.023zA) 

= _^™ fcA 7 2 + A7 4 (^ fc arccos I + ^A) ee iS^ T , (42) 

where the expression in the last line is recognized as the standard form of the action 
for Euclidean dynamical triangulations in four dimensions. 
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5 The transfer matrix 



We have now all the necessary prerequisites to study the full partition function or 
amplitude for pure gravity in the dynamically triangulated model, 

£ cm e ' s<T> ' (43) 

triangulations T 

where the sum is taken over inequivalent Lorentzian triangulations, corresponding 
to inequivalent discretized space-time geometries. We have chosen the measure 
factor to have the standard form in dynamical triangulations, m(T) = 1/C(T). 
The combinatorial weight of each geometry is therefore not 1, but the inverse of 
the order C(T) of the symmetry (or automorphism) group of the triangulation T. 
One may think of this factor as the remnant of the division by the volume of the 
diffeomorphism group Diff(M) that would occur in a formal gauge-fixed continuum 
expression for Z. Its effect is to suppress geometries possessing special symmetries. 
This is analogous to what happens in the continuum where the diffeomorphism 
orbits through metrics with special isometries are smaller than the "typical" orbits, 
and are therefore of measure zero in the quotient space Metrics/Diff(M). 

Again, we expect the final result to be independent of the details of the dynamical 
and combinatorial weights in (|i"3D. The sparseness of critical points and universality 
usually ensure that the same continuum limit (if existent) is obtained for a wide 
variety of initial regularized models. 

The most natural boundary conditions for Z in our discretized Lorentzian model 
are given by specifying spatial discretized geometries g at given initial and finite 
proper times t. The slices at constant t are by construction space-like and of fixed 
topology, and are given by (unlabelled) Euclidean triangulations in terms of {d — 1)- 
dimensional building blocks (equilateral Euclidean triangles in d = 3 and equilateral 
Euclidean tetrahedra in d — 4) . 

Given an initial and a final geometry gi and g2, we may think of the corre- 
sponding amplitude Z as the matrix element of the quantum propagator of the 
system, evaluated between two states \gi) and \g 2 ). Since we regard distinct spatial 
triangulations gi as physically inequivalent, a natural scalar product is given by 

(0ita> = £^<W ^20{g)\g){g\ = l (44) 

In line with our previous reasoning, we have included a symmetry factor C(g) for 
the spatial triangulations. 

In the regularized context, it is natural to have a cutoff on the allowed size 
of the spatial slices so that their volume is < N. The spatial discrete volume 
Vol(g) is simply the number of (d — l)-simplices in a slice of constant integer t. We 
define the finite-dimensional Hilbert space as the space spanned by the vectors 
{\g), N min < Vol(g) < N}, endowed with the scalar product (fP]). The lower bound 
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iV m i n is the minimal size of a spatial triangulation of the given topology satisfying 
the simplicial manifold conditions. It is well-known that the number of states in 
H( N ) is exponentially bounded as a function of N [E4[ . 



For given volume-cutoff N, we can now associate with each time-step At = 1 a 
transfer matrix Tpj describing the evolution of the system from t to t+1, with matrix 
elements 

(g 2 \f N (a)\g 1 )=G a (g 1 ,g 2 ;l) = J] _L e^-CO ^ 

The sum is taken over all distinct interpolating <i-dimensional triangulations g\ — > g 2 
of the "sandwich" with boundary geometries g\ and g 2 , contributing AS to the 
action, according to (|35|) , (|40|). The propagator Gj^(gi, g 2 ;t) for arbitrary time 
intervals t is obtained by iterating the transfer matrix t times, 

G N (g 1 ,g 2 ;t) = (g 2 \f t N \g 1 ), (46) 

and satisfies the semigroup property 

G N (g u g 2 ; h + t 2 ) = ^ C(g) G N (g u g; t x ) G N (g, g 2 ; t 2 ), (47) 

g 

where the sum is over all spatial geometries (of bounded volume) at some interme- 
diate time. Because of the appearance of different symmetry factors in ( |45|) and 
( f47j ) it is at first not obvious why the composition property ([47]) should hold. In 
order to understand that it does, one has to realize that by continuity and by virtue 
of the manifold property there are no non-trivial automorphisms of a sandwich 
T = gi — > g 2 that leave its boundary invariant. It immediately follows that the 
(finite) automorphism group of T must be a subgroup of the automorphism groups 
of both of the boundaries, and that therefore C(T) must be a divisor of both C(g\) 
and C(g 2 ). It is then straightforward to verify that the factor C(g) appearing in the 
composition law ([17]) ensures that the resulting geometries appear exactly with the 
symmetry factor they should have according to ([45|). 

Our construction of the propagator obeying ([IT]) parallels the two-dimensional 
treatment of @. Also in higher dimensions one may use the trick of introducing 
"markings" on the boundary geometries to absorb some of the symmetry factors C. 
However, note that unlike in d = 2 the number of possible locations for the marking 
does not in general coincide with the order C(g) of the boundary's symmetry group. 



6 Properties of the transfer matrix 

In this section we will establish some desirable properties of the transfer matrix 
Tjy(a), defined in ([45]), which will enable us to define a self-adjoint Hamiltonian op- 
erator. Since the name "transfer matrix" in statistical mechanical models is usually 
reserved for the Euclidean expression T = e~ ah , where a denotes the lattice spacing 
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in time-direction, we will from now on work with the Wick-rotated matrix Tjv(— 1 



A set of sufficient conditions on Tjy guaranteeing the existence of a well-defined 
quantum Hamiltonian h are that the transfer matrix should satisfy 

(a) symmetry, that is, T N = Tjy. This is the same as self-adjointness since the 
Hilbert space which Tjv acts on is finite-dimensional. It is necessary if 
the Hamiltonian is to have real eigenvalues. 

(b) Strict positivity is required in addition to (a), that is, all eigenvalues must be 
greater than zero; otherwise, tiN = — cT 1 logT/v does not exist. 

(c) Boundedness, that is, in addition to (a) and (b), T N should be bounded above 
to ensure that the eigenvalues of the Hamiltonian are bounded below, thus 
defining a stable physical system. 

Establishing (a)-(c) suffices to show that our discretized systems are well-defined as 
regularized statistical models. This of course does not imply that they will possess 
interesting continuum limits and that these properties will necessarily persist in the 
limit (as would be desirable). On the other hand, it is difficult to imagine how 
the continuum limit could have such properties unless also the limiting sequence of 
regularized models did. 

All of the above properties are indeed satisfied for the Lorentzian model in d = 2 
|7j, where moreover the quantum Hamiltonian and its complete spectrum in the 



continuum limit are known explicitly |25|, 18| . Note that self-adjointness of the 



continuum Hamiltonian H implies a unitary time evolution operator e~ tHT if the 
continuum proper time T is analytically continued. In d > 2 we are able to prove 
a slightly weaker statement than the above, namely, we can verify (a)-(c) for the 
two-step transfer matrix T^.[] This is still sufficient to guarantee the existence of a 
well-defined Hamiltonian. 

One verifies the symmetry of the transfer matrix by inspection of the explicit 
form of the matrix elements (0). The "sandwich actions" AS a as functions of the 
boundary geometries gx, g<i in three and four dimensions can be read off from (]35|) 
and (f4Up. To make the symmetry explicit, one may simply rewrite these actions as 
separate functions of the simplicial building blocks and their mirror images under 
time-reflection (in case they are different). Likewise, the symmetry factor C(T) 
and the counting of interpolating geometries in the sum over T are invariant under 
exchange of the in- and out-states, \gi) and | ^72 ) - 



Next, we will discuss the reflection (or Osterwalder-Schrader) positivity [26. 27 



of our model, with respect to reflection at planes of constant integer and half-integer 



time (see also |28| and references therein). These notions can be defined in a straight- 
forward way in the Lorentzian model because it possesses a distinguished notion of 
(discrete proper) time. Reflection positivity implies the positivity of the transfer 
matrix, TV > 0. 



7 This corrects a claim in [p5| where we stated that (b) holds for the one-step transfer matrix. 
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"Site reflection" denotes the reflection 6 S with respect to a spatial hypersurface 
of constant integer-t (containing "sites", i.e. vertices), for example, t = 0, where it 
takes the form 

6 S : t -> -t. (48) 

Let us accordingly split any triangulation T along this hypersurface, so that T~ is 
the triangulation with t < and T + the one with t > 0, and T _ flT + = g(t = 0) = (70, 
where go denotes a spatial triangulation at t = 0. Consider now functions F that 
depend only on T + (that is, on all the connectivity data specifying T + uniquely, in 
some parametrization of our choice). Site-reflection positivity means the positivity 
of the Euclidean expectation value 

({6 8 F)F) > 0, (49) 

for all such functions F . The action of 9 S on a function F(T + ) is defined by anti- 
linearity, (6 S F)(T~) := F{6 S (T~)). By virtue of the composition property ((4?D, we 
can write 



((dsF)F) = ^E^)(W 



e ~S(T) 



-S(T+) 



go t- T+ 

r-ct=o)= 9o T+(i=o)= 90 

= r^^ojVW^O. (50) 
90 

The equality in going to the last line holds because both the action and the symmetry 
factor C depend on "outgoing" and "incoming" data in the same way (for example, 
on (m,n)-simplices in the same way as on (n,m)-simplices). Note that the standard 
procedure of extracting a scalar product and a Hilbert space from ([50]) is consistent 
with our earlier definition (0). One associates (in a many-to-one fashion) functions 
F(T + ) with elements of a Hilbert space at fixed time t with scalar product (•, •), 
where ^g) = ((@sF)G) p6| . A set of representative functions which reproduces 
the states and orthogonality relations ( fP]) is given by 



+)H Vv^), n<* = 0) = 9 

otherwise, 

as can be verified by explicitly computing the expectation values ((6 s F g i)F g ii) . 

We have therefore proved site-reflection positivity of our model. This is already 
enough to construct a Hamiltonian from the square of the transfer matrix T/y (see 
eq. (^) below), since it implies the positivity of Tfj E8|]. 

Proving in addition link-reflection positivity (which would imply positivity of the 
"elementary" transfer matrix, and not only of Tjy-) turns out to be more involved. A 
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"link reflection" is the reflection 9 t at a hypersurface of constant half-integer-t, for 
example, t = 1/2, 

0, : t^l-t. (52) 
To show link-reflection positivity in our model we would need to demonstrate that 

((OiF)F) > 0, (53) 

where F is now any function that depends only on the part T + of the triangulation 
T at times later or equal to 1. We can again write down the expectation value, 



((9iF)F) = Z- 1 J2Y, C ^o)C(9i)G N (g ,g 1 ;l) 



x 

90 91 



y- (0 s F)(T ) s(T - } F(T+) 5(T+) 

2L, c(t-) ^ c(T+) e • ^ 

T-(t=0)=g T+(t=l)= sl 

In order to show that this is positive, one should try to rewrite the right-hand side as 
a sum of positive terms. The proof of this is straightforward in d — 2 (see App. 2 for 
details), but it is considerably more difficult to understand what happens in higher 
dimensions. The reason for this is the non-trivial way in which the various types 
of simplicial building blocks fit together in between slices of constant integer-t. In 
a way, this is a desirable situation since it means that there is a more complicated 
"interaction" among the simplices. It is perfectly possible that T N itself is not 
positive for d = 3,4. This may depend both on the values of the bare couplings and 
on the detailed choices we have made as part of our discretization. (By contrast, 
it is clear from our proof of site-reflection positivity that this property is largely 
independent of the choice of building blocks.) 

Nevertheless, as already mentioned above, site-reflection positivity is perfectly 
sufficient for the construction of a well-defined Hamiltonian. So far, we have only 
shown that the eigenvalues of the squared transfer matrix are positive. In order for 
the Hamiltonian 

h' N := ~]ogTfr (55) 
to exist, we must achieve strict positivity. We do not expect that the Hilbert space 

H (N) 

contains any zero-eigenvectors of T/v since this would entail a "hidden" sym- 
metry of the discretized theory. It is straightforward to see that none of the basis 
vectors \g) can be zero-eigenvectors. However, we cannot in principle exclude "acci- 
dental" zero-eigenvectors of the form of linear combinations Y^i a i\9i)- m case sucn 
vectors exist, we will simply adopt the standard procedure of defining our physical 
Hilbert space as the quotient space H^p = / 'Af( N \ where Af( N ^ denotes the 
span of all zero-eigenvectors. 

Lastly, the boundedness of the transfer matrix (and therefore also of Tjj) for 
finite spatial volumes N follows from the fact that (i) there is only a finite number 
of eigenvalues since the Hilbert space is finite-dimensional, and (ii) each matrix 
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element (<72|Tjv|<7i) has a finite value, because it has the form of a finite sum of terms 
e~ s /C(T). (Note that this need not be true in general if we abandoned the simplicial 
manifold restrictions, because then the number of interpolating geometries for given, 
fixed boundary geometries would not necessarily be finite.) 



7 Monte Carlo moves 

In numerical simulations of the discrete Lorentzian model, as well as for some the- 
oretical considerations, one needs so-called "moves", that is, a set of basic (and 
usually local) manipulations of a triangulation, resulting in a different triangulation 
within the class under consideration. The relevant class of geometries for our pur- 
poses is that of Lorentzian (three- or four-dimensional) triangulations with a fixed 
number t of time-slices, with finite space-time volume, and obeying simplicial man- 
ifold constraints. The set of basic moves should be ergodic, that is, one should be 
able to get from any triangulation to any other by a repeated application of moves 
from the set. 

The Monte Carlo moves used in Euclidean dynamical triangulations |2^, [5] are 
not directly applicable since in general they do not respect the sliced structure of 
the Lorentzian discrete geometries. Our strategy for constructing suitable sets of 
moves is to first select moves that are ergodic within the spatial slices t = const 
(this is clearly a necessary condition for ergodicity in the full triangulations), and 
then supplement them by moves that act within the sandwiches At = 1. In d = 3, 
this results in a set of five basic moves which has already been used in numerical 



simulations [|1^, |30|, and which we believe is ergodic. For completeness, we will 
describe them again below. We will also propose an analogous set of moves in 
four dimensions, which are somewhat more difficult to visualize. In manipulating 
triangulations, especially in four dimensions, a useful mental aid are the numbers of 
simplicial building blocks of type A contained in a given building block of type B of 
equal or larger dimension (Table 1). 

As usual, if the moves are used as part of a Monte Carlo updating algorithm, 
they will be rejected in instances when the resulting triangulation would violate the 
simplicial manifold constraints. The latter are restrictions on the possible gluings 
(the pair- wise identifications of sub-simplices of dimension d—1) to make sure that 
the boundary of the neighbourhood of each vertex has the topology of a 5' d ~ 1 -sphere 
|[23|| . In dimension three, "forbidden" configurations are conveniently characterized 
in terms of the intersection pattern of the triangulation at half-integer times (see 
|T5| |, Appendix 2). This is given by a two-dimensional tesselation in terms of squares 
and triangles or, alternatively, its dual graph with three- and four-valent crossings. 
Also in four dimensions, one may characterize sandwich geometries by the three- 
dimensional geometry obtained when cutting them at height 1/2. It is formed by 
two types of three-dimensional building blocks, a tetrahedron and a "thickened" 
triangle. 
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Table 1: Numbers of simplicial building blocks contained in simplices of equal or 
higher dimension. 

7.1 Moves in three dimensions 

As already described in |TB[, there are five basic moves (counting inverse moves as 
separate). They map one 3d Lorentzian triangulation into another, while preserving 
the constant-time slice structure, as well as the total proper time t. We label the 
moves by how they affect the number of simplices of top-dimension, i.e. d = 3. A 
(m, n)-move is one that operates on a local sub-complex of m tetrahedra and replaces 
it by a different one with n tetrahedra. The tetrahedra themselves are characterized 
by 4-tuples of vertex labels. Throughout, we will not distinguish moves that are 
mirror images of each other under time reflection. In detail, the moves are 

(2,6): this move operates on a pair of a (1,3)- and a (3,l)-tetrahedron (with vertex 
labels 1345 and 2345) sharing a spatial triangle (with vertex labels 345). A 
vertex (with label 6) is then inserted at the centre of the triangle and connected 
by new edges to the vertices 1, 2, 3, 4 and 5. The final configuration consists 
of six tetrahedra, three below and three above the spatial slice containing the 
triangles (Fig. [|). This operation may be encoded by writing 

1 345 + 2345 -> 1346 + 2346 + 1356 + 2356 + 1456 + 2 456 , (56) 

where the underlines indicate simplices shared by several tetrahedra (not nec- 
essarily all, in order not to clutter the notation). The inverse move (6,2) 
corresponds to a reversal of the arrow in fl5"6"|). Obviously, it can only be per- 
formed if the triangulation contains a suitable sub-complex of six tetrahedra. 
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Figure 3: The (2,6)-move in three dimensions. 




Figure 4: The (4,4)-move in three dimensions. 

(4,4): This move can be performed on a sub-complex of two (1,3)- and two (3,1)- 
tetrahedra forming a "diamond" (see Fig. [|), with one neighbouring pair each 
above and below a spatial slice. The move is then 

1 235 + 2356 + 1345 + 3456 -> 1234 + 2346 + 1245 + 2456. (57) 

From the point of view of the spatial "square" (double triangle) 2345, the move 
(|57j) corresponds to a flip of its diagonal. It is accompanied by a corresponding 
reassignment of the tetrahedra constituting the diamond. The (2,6)- and (6,2)- 
moves, together with the (4,4)-move (which is its own inverse) induce moves 
within the spatial slices that are known to be ergodic for two-dimensional 
triangulations. 

(2,3): The last move, together with its inverse, affects the sandwich geometry without 
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Figure 5: The (2,3)-move in three dimensions. 



changing the spatial slices at integer-t. It is performed on a pair of a (3,1)- 
and a (2,2)-tetrahedron which share a triangle 345 in common (see Fig. |5|), 
and consists in substituting this triangle by the one- dimensional edge 12 dual 
to it, 

1345 + 2345 -^ 1234+ 1235 + 1245. (58) 

The resulting configuration consists of one (3,1)- and two (2,2)-tetrahedra, 
sharing the link 12. Again, there is an obvious inverse. 

7.2 Moves in four dimensions 

If we distinguish between the space- and time-like character of all of the four- 
dimensional moves, there is a total of ten moves (including inverses). We will again 
characterize simplices in terms of their vertex labels. The first two types of moves, 
(2,8) and (4,6), reproduce a set of ergodic moves in three dimensions when restricted 
to spatial slices. We will describe each of the moves in turn. 

(2,8): The initial configuration for this move is a pair of a (1,4)- and a (4,l)-simplex, 
sharing a purely space-like tetrahedron 3456 (Fig. |]). The move consists in 
inserting an additional vertex 7 at the centre of this tetrahedron and subdiving 
the entire sub-complex so as to obtain eight four-simplices, 

13456+23456 -> 13457+23457+13467+23467+13567+23567+14567+24567, 

(59) 

with an obvious inverse. 

(4,6): In this configuration, we start from a pair (2345,3456) of spatial tetrahedra 
sharing a common triangle 345, which are connected to a vertex 1 at time 
t—1 and another vertex 7 at time t+1, forming together a sub-complex of size 
six (Fig. |7|). The move consists in swapping the triangle 345 with its (spa- 
tially) dual edge 26, during which the two spatial tetrahedra are substituted 
by three, and the geometry above and below the spatial slice at time t changed 
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Figure 6: The (2, 8) -move in four dimensions. 

accordingly. In our by now familiar notation, this amounts to 

12345+23457+13456+34567 -> 12346+23467+12356+23567+12456+24567, 

(60) 

with the arrow reversed for the corresponding inverse move. 

This type of move comes in two varieties. Its general structure is as follows: 
the initial configuration is a pair of 4-simplices with a common tetrahedron. 
During the move, this tetrahedron is "deleted" and substituted by its (in a 
four-dimensional sense) dual edge. The end result is a sub-complex consisting 
of four 4-simplices. From a Lorentzian point of view, there are two situations 
where the application of this move does not interfere with the slice-structure 
or the manifold constraints. In the first one, a (4,1)- and a (3,2)-tetrahedron 
from the same sandwich share a (3,l)-tetrahedron 3456 (Fig. |8|). The dual edge 
12 is time-like and shared in the final configuration by one (4,1)- and three 
(3,2)-simplices. The second possibility is that of two (3,2)-simplices sharing 
a (2,2)-tetrahedron. One of the (3,2)-simplices is "upside-down", such that 
the entire sub-complex has spatial triangles in both the slices at t and at t+1 
(134 and 256, see Fig. |). After the move, the total number of (3,2)-simplices 
has again increased by two. (Note that the changes of the bulk variables are 
encoded in the /-vectors of Sec. ||.) Both types of moves are described by the 
relation 

1 3456 + 2 3456 -» 12345 + 12346 + 12456 + 12356, (61) 
and their inverses by the converse relation. 

The initial sub-complex in this type of move is made up of three 4-simplices 
which share a triangle in common. In the course of the move, this triangle is 
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Figure 7: The (4,6)-move in four dimensions. 




Figure 8: The (2,4)-move in four dimensions, first version. 

"deleted" and substituted by its dual (in a four- dimensional sense), which is 
again a triangle. It is straightforward to verify that this move can only occur in 
Lorentzian gravity if both of the triangles involved are time-like. Again, there 
are two allowed variations of the move. In the first one, both the initial and 
final sub-complex consist of one (4,1)- and two (3,2)-simplices, and the space- 
like edge of each of the triangles 123 and 456 lies in the same slice t = const 
(Fig. [TI|). The 4-simplices are rearranged according to 

12456 + 13456 + 23456 -> 12345 + 12346 + 12356. (62) 

The other initial configuration for which flB2|) can be performed involves only 
(3,2)-simplices of both orientations, as illustrated in Fig. [11]. The "swapped" 
triangle 123 now has its space-like edge in the opposite spatial slice from that 
of the original triangle 456. As in the (4,4)-move in three dimensions, this 
type of move is its own inverse. 
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Figure 9: The (2,4)-move in four dimensions, second version. 




5 5 



Figure 10: The (3,3)-move in four dimensions, first version. 

8 Kinematic bounds 

Without solving the combinatorics of the Lorentzian gravity models explicitly, some 
qualitative insights can be gained by studying the behaviour of their partition func- 
tions for either very small or very large values of the inverse Newton's constant k. By 
comparing our results with those of Euclidean dynamical triangulations in 3 and 4d 
3TJ, 0, we will see that significant differences arise already at this kinematic level. 



Since the partition function for Lorentzian gravity is evaluated after Wick-rotating, 
we will in the following set a = — 1, and work in the Euclidean sector of our models. 
We will perform the analysis for d = 3 and d = 4 in turn. 

8.1 Analysis in three dimensions 

For ease of comparison with the Euclidean results, we write the Regge action in 
standard form (c.f. ([57]) ), namely, as 

4dt = hN 3 - hN h (63) 
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Figure 11: The (3,3)-move in four dimensions, second version. 



with the couplings 

ki = 2nk, k 3 = 6 k arccos — h A — ■=. (64) 

3 6^ 

Lorentzian triangulations T are conveniently characterized by their /-vectors f(T), 
which keep track of the numbers (|T9| ) of simplicial building blocks. In three dimen- 
sions, we have 

/( 3 ) = (JVq, N**, iV x TL , iVf , iV 2 TL , NjV\ Nt 2) ). (65) 
Rewriting the action as 

Sg T = N 3 (k 3 - k^) =: N 3 (k 3 - (66) 
iV 3 

we see that for large coupling k\ and at fixed space-time volume, configurations with 
a large value of £ will be energetically preferred. In the Euclidean case, it can be 
shown that in the thermodynamic limit ^3 — ► 00, the ratio £ obeys the kinematic 
bounds 

4 

1 < £ < -, Euclidean 3d triangulations. (67) 
3 

The upper bound is saturated for so-called branched-polymer (or "stacked-sphere" ) 
configurations. Geometrically degenerate configurations of this type dominate the 
partition function in the continuum limit and seem to be responsible for the failure 
of the Euclidean model to lead to an interesting continuum theory, both in three 
and four dimensions. 

Obviously, the geometries of the Lorentzian model must satisfy flB7p, which is 
valid for all 3d simplicial manifolds. However, it turns out that the Lorentzian 
dynamical triangulations satisfy a more stringent upper bound. Using identities 
from Sec. |3|, we can write 

N 3 = Nj?' 2) + NjfV = Af ' 2) + 4iV - A X t > + 4iVo _ 4 x t. (68) 
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The minimal A^ 2 ' 2 -* is given by N^ 2 ^ n = c m i n t, where c m i n > is an integer counting 
the minimal possible number of (2,2)-tetrahedra per time-step At — 1. Using further 
identities, fl68| ) therefore becomes 

4iV x < 5N 3 - (c min - A X )t. (69) 

From this we deduce that in the large- volume limit (where t/N 3 — > 0) the upper 
bound for £ is lowered to 5/4, so that the Euclidean relation ( |6^) is substituted by 

5 

1 < £ < — , Lorentzian 3d triangulations. (70) 

We conclude that the branched-polymer regime of Euclidean gravity cannot be real- 
ized in the Lorentzian case. The upper bound in Q7DD is saturated when the number 
of (2,2)-simplices is minimal. 

There is an alternative way of deriving ([701), which also sheds some light on 
the nature of the lower bound. Here, one asks whether the bounds on £ can be 
approached by a successive application of the Monte Carlo moves described in Sec. 
0. In three dimensions, the effect of the Monte Carlo moves on the /-vector is 

A (2 ,6)/ (3) = (1,3,2,2,6,4,0), (71) 
A (4 ,4)/ (3) = (0,0,0,0,0,0,0), (72) 
A (2 ,3)/ (3) = (0,0,1,0,2,0,1), (73) 

with obvious expressions for the inverse moves. Let us now choose some initial 
configuration /q , for example, a t-translation-invariant geometry of minimal volume 

(3) -» -» 

of the form /q = c i, where c is a constant seven- vector characterizing a sandwich 
geometry. Any triangulation that is reached from the initial one by a sequence of 
allowed Monte Carlo moves has an /-vector 

/( 3 ) = (at + x, c 2 t + 3x, c 3 t + 2x + y, c 4 t + 2x, c 5 t + Qx + 2y, c 6 t + Ax, c 7 t + y). (74) 

The integers x = n(2,6) — n (6,2) and y = n^^—n^) count the "excess" of applications 
of a move over its inverse. From ([74D we can read off the value of £ as 



. _ Ni _ (c 2 + c 3 )t + 5x + y 
N 3 (c 6 + c 7 )t + Ax + y 

It is obvious that for large triangulations this ratio is maximized by applying only 
(2,6)-moves. This is always possible, since this Monte Carlo move is a subdivision 
move that can always be performed. In the infinite-volume limit (this also implies 
x 3> t), £ approaches 5/4, in accordance with our earlier result. Alternatively, 
a triangulation can be enlarged by applying (2,3)-moves, which according to (|73|) 
increase the number of (2,2)-tetrahedra. It is not immediately clear that we can 
saturate the lower bound in (^) by taking y — > oo since y is not entirely inde- 
pendent of x (the applicability of the (2,3)-move is related to the availability of 
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(3,l)-tetrahedra, which in turn is governed by x). Nevertheless, in the large- volume 
limit such configurations do actually exist, as we will demonstrate below. 

Another type of degenerate geometry that was observed to dominate the state 
sums of Euclidean dynamical triangulations at small k []33], |34| is characterized by 
an inequality which is non-linear in the bulk variables, namely, 

Nt < ^N (N - 1). (76) 



A triangulation saturating ( fT6|) is called "2-neighbourly" |32| since any of its vertices 



is connected by an edge to any other vertex. The resulting space is of maximal 
Hausdorff dimension, because "one can get anywhere in a single step" . Such singular 
structures have indeed been observed in the so-called crumpled phase of Euclidean 
quantum gravity in three and four dimensions. 

Again, it is impossible to come anywhere near this phase in the continuum limit 
of the Lorentzian model. This has to do with the fact that by construction of the 
Lorentzian geometries, two vertices can only be connected by a link if they lie in 
the same or in successive spatial slices. Instead of (|T6]), we therefore have separate 
relations for the numbers Nf L and N± L of space- and time-like edges, 



t-i t-i 



N? L = ^(3iV (r) - 6) = 3N - 6t, N? L < £ N (t)N (t + 1). (77) 



Assuming canonical scaling, the right-hand side of the inequality ( [76]) behaves like 
(length) 6 , whereas the second relation in Q77D scales only like (length) 5 , proving our 
earlier statement. 

The inequality in ( fT7|) can indeed be saturated by certain regular Lorentzian 
geometries in the large- volume limit, as one can show by explicit construction. For 
ease of visualization, let us assume that the two-dimensional slices have the topology 
of two-tori T 2 . Consider a sandwich where the initial two-geometry at time r is a 
"strip" of discrete length li and width c\ = 2, consisting of 1l\ + 4 triangles^, with 
opposite sides identified to give a torus. Similarly, at time r + 1, take a strip of 
length I2 and width C2 = 2 with 2/2 + 4 triangles. Choose the interpolating geometry 
between the two tori such that the (2,2)-tetrahedra when viewed at time r+1/2 form 



two rectangular blocks, of size l-J, 2 and of size Cic 2 (Fig. [12]) . Consider now a sequence 
of such configurations where both of the li become large, but the integers q remain 
unchanged. Because of N q (t) = Nf h (r)/3 = h + 2 and JV (r + l) = JV? L (r + l)/3 = 
I2 + 2, we find that in the large- volume limit the bulk variables associated with the 
sandwich are given by Aq « N? L « hl 2 , and by N 3 = A^ 2 SL (r) + iYf L (r+ 1) + A^ 3 (2 ' 2) w 
^1^2- For the vertex numbers, we have N (t)N (t + 1) m l x l 2 . This verifies our earlier 
claims that both the inequality ([77|) and the lower bound in flTOj ) are saturated by this 

8 The +4 comes from the fact that the minimal length of the intermediate layer of the two-torus 
must be I = 2, to give a simplicial 2d manifold under gluing. 
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Figure 12: One of a sequence of toroidal sandwich configurations at r + 1/2 which 
saturate the inequality in flTT|) in the large-volume limit (opposite sides to be identi- 
fied). Light links correspond to triangle edges at time r, dark links to triangle edges 
at time r + 1. 



class of configurations. Of course, we achieved this by a "fine-tuning" of geometry, 
and the configurations thus obtained are highly non-generic. 

The above considerations show that the phase structure of three-dimensional 
Lorentzian gravity must be quite different from that of the Euclidean theory since 
already at a kinematic level, the extreme branched-polymer and crumpled configu- 
rations cannot be realized. Obviously it does not exclude the possibility that milder 
pathologies prevail in the continuum limit, but it is clearly a step in the right direc- 
tion. Hopes that the Lorentzian model is therefore better-behaved than its Euclidean 
counterpart have been boosted by recent numerical investigations in three dimen- 
1TEL |3U| . There it was found that for a large range of the gravitational coupling 



sions 



constant k\, the partition function is dominated by a phase of extended geometry 
(and the parameter £ stays well away from its extremal bounds ([70"D), whose global 
scaling properties are those of a three-dimensional object. Compared with the Eu- 
clidean case, this is a completely new and very promising phenomenon. It reiterates 
the finding in d — 2 J7| that the causal structure of the Lorentzian space-times acts 
effectively as a regulator for the quantum geometry. 

8.2 Analysis in four dimensions 

The analysis of the kinematic bounds in four dimensions is completely analogous to 
that of the previous subsection, and we will refrain from spelling out all the details. 
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We start from the action in the form (c.f. (0)) 

Sg> T = k 4 N, - k 2 N 2 = N 4 (k 4 - fc3) =: N 4 (k 4 - feaO, (78) 
with the coupling constants 

„ = ^ rt , fa = M, i + A ±§. (79) 

2 z 4 96 



In (|78f) we have already anticipated that a convenient order parameter in four di- 
mensions is the quotient £' = N 2 /N i . The /-vector has now length 10, and is given 
by 

/(*) = (JVo, iVf L , jVf-, iVf, iV 2 TL , iVf, jvf > x \ ivf ' 2) , jvi 4ll) , i\f ' 2) ). (80) 

The discussion parallels closely the three-dimensional case, with the relevant Eu- 
clidean inequalities []3l], [32| 



5 

2 < £' < — , Euclidean 4d triangulations. (81) 

Let us derive their Lorentzian analogues by analyzing the effect of the Monte Carlo 
moves of Sec. [7|on the /-vector flS0|), which is 

A (2i8) /( 4 ) = (1,4,2,6,8,3,12,0,6,0), (82) 

A(4, 6 )/ {4) = (0,1,0,2,2,1,4,0,2,0), (83) 

A (2 ,4)/ (4) = (0,0,l,0,4,0,2,3,0,2), (84) 

A (3 ,3)/ (4) = (0,0,0,0,0,0,0,0,0,0). (85) 

Starting from a minimal initial configuration f{f \ and defining x = ri(2,8) — ^(8,2), 
V = n (4,6) — ra (6,4) and z = ri(2,4) — ra (4,2), and dropping terms proportional to t in the 
large- volume limit, one obtains 

£' = ^ — > 7x + 2(^ + 2;) 
A^ 4 3x + y + 2; 

Again, £' is maximized by applying the subdivision move (2,8) only. Note also that - 
like in three dimensions - this corresponds to creating stacked-sphere configurations 
within the spatial slices of constant integer-t. (We are using the expression "stacked 
spheres" in a loose sense, and only require that the triangulations look like stacked 
spheres almost everywhere.) For the purposes of relation flSBp, the moves (4,6) and 
(2,4) play an equivalent role, although they generate different types of four-simplices. 
It is clear that a lower bound for the ratio is given by 2, although again it is not 
immediately clear that this value can actually be attained by any triangulation in the 
thermodynamic limit. Like in 3d, we can conclude that the range for the parameter 
£' in Lorentzian gravity is reduced, such that 

7 

2 < £ < — , Lorentzian 4d triangulations. (87) 
3 
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We leave it as an exercise to the reader to prove this relation directly from (p7\) and 
the identities in Sec. 3, as well as to construct explicit triangulations that saturate the 
lower bound. Again, our analysis implies the absence of the polymeric geometries 
found in the Euclidean theory. (Note that a polymeric behaviour is in principle 
allowed both in three- and four- dimensional Lorentzian triangulations, as long as it 
takes place in the spatial directions only. In fact, this corresponds precisely to the 
limit x — > oo in both (|75|) and (|86|).) Likewise, there is no crumpled phase. The 
argument of d = 3 carries over to four dimensions, with the only difference that the 
right-hand side of (|76|) scales like (length) 8 , whereas the corresponding relation in 
(|77|) scales like (length) 7 . 



9 Summary and outlook 

In this paper, we have presented a comprehensive analysis of a regularized path 
integral for gravity in terms of Lorentzian dynamical triangulations. We have shown 
that for a finite space-time volume N the model is well-defined as a statistical system 
of geometry, with a finite state sum. 

Having established these properties, our main interest is the behaviour of the 
system in a suitable continuum limit. For this, we must take the lattice volume 
N — > oo and the invariant lattice spacing a — > 0, and allow for a renormalization of 
both the bare coupling constants and the propagator itself. Although the discretized 
set-up of our model looks rather similar in different dimensions, we expect the details 
of the renormalization and the resulting physics to depend strongly on d. 

The continuum theory in d = 2 is already explicitly known, since the Lorentzian 
model has been solved exactly [^, ||. It is obtained by simply taking iV — > oo or, 
equivalently, tuning the bare cosmological constant A to its critical value. The theory 
describes a one-dimensional spatial "quantum universe" whose volume fluctuates in 
time. Obviously the physical significance of this model is somewhat obscure, since 
there is no classical theory of 2d Einstein gravity, and therefore no well-defined 
classical limit. Nevertheless, it provides a beautiful example of the inequivalence of 
path integrals of Euclidean and Lorentzian geometries. What is more, we understand 



in a very explicit manner |7|, |35| that the origin of the difference lies in the presence 
and absence of "baby universes" branching off in the time- direct ion (which in the 
Lorentzian case are incompatible with the "micro- causality" of the path-integral 
histories) . 

"Solving the path integral" in our formulation amounts to solving the combina- 
torics of the Lorentzian triangulations in a given dimension and for given boundary 
data, at least asymptotically. Although this counting problem is mathematically 
well-defined, it is not simple when d > 3. In three space-time dimensions, where 
the problem is presently under investigation, an interesting relation with an already 
known matrix model with ABAB-interaction has been uncovered [|l9fl . 

A potentially powerful feature of our models is the fact that they can also be 
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studied numerically. Using Monte Carlo methods, it was shown that 2d Lorentzian 
gravity coupled to matter has no barrier at c = 1 , and that pure 3d Lorentzian 
gravity possesses a ground state of extended geometry, with the global scaling prop- 
erties of a three-dimensional space-time universe [|l6i |30|| . Both of these findings are 
highly significant, because they show that these models lie in different universality 
classes from their Euclidean counterparts. 

Let us recapitulate how this difference with the purely Euclidean formulation 
came about. Our physical premise was to take the space of Lorentzian space-times 
as our starting point. This enabled us to impose certain causality constraints on 
the individual histories contributing to the path integral. As a consequence, like in 
classical relativity, "time"- and "space" -directions are not interchangeable, a feature 
that persists in the continuum quantum theory. In order for this philosophy to lead 
anywhere (i.e. to yield convergent mathematical expressions) we crucially needed a 
Wick rotation on the space of Lorentzian geometries. The fact that such a map exists 
in our formulation is by no means trivial. For example, if a similar construction were 
attempted in Regge calculus, one would immediately run into trouble with triangular 
inequalities. 

Being familiar with the pitfalls of Euclidean gravitational path integrals, one 
may now wonder whether our path integrals after the Wick rotation to Euclidean 
signature do not also suffer from pathologies, in the limit as N — > oo. In particular, 
one may worry about the presence of a conformal divergence in d > 3 associated 
with the unboundedness of the gravitational action. Remarkably, although such a 
divergence is potentially present in our formulation, since configurations with large 
negative action exist, these are entropically suppressed in the continuum limit, at 



least in d = 3 It has been demonstrated in |37[ how this phenomenon can be 



understood from the point of view of the continuum gravitational path integral: a 
careful analysis of various Faddeev-Popov determinants reveals how the resulting 
path-integral measure can lead to a non-perturbative cancellation of the conformal 
divergence^ 

To summarize, we have constructed a new and powerful calculational tool for 
quantum gravity that in the best of all worlds will lead to a non-perturbative defini- 
tion of an interacting quantum theory of gravity in four dimensions. An inclusion of 
other matter degrees of freedom is in principle straightforward. At the discretized 
level, our models are as well-defined as one could hope for. An investigation in 
dimension 3, both analytically and numerically, is well under way. It remains a ma- 
jor challenge to set up sufficiently large numerical simulations in four dimensions. 
Clearly, because of the anisotropy between time and spatial directions, larger lattice 
sizes than in the Euclidean dynamically triangulated models will be needed. On the 
analytic front, tackling the full combinatorics in four dimensions will be difficult, as 
one would expect from any fully-fledged formulation of quantum gravity. In order to 
approach the problem, we are currently studying both the three-dimensional theory 

9 Note that the use of proper-time gauge in J3?|] made it possible to define a non-perturbative 
Wick rotation analogous to the one used in the present work. 
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and its relation to other formulations of 3d quantum gravity, as well as gravitational 
models with extra symmetries. 
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Figure 13: Positive (a) and negative (b) Euclidean deficit angles S. 

Appendix 1: Lorentzian angles 

For the sake of completeness, we will describe some properties of Lorentzian an- 
gles (or "boosts") which appear in the Regge scalar curvature as rotations about 
space-like bones (that is, space-like links in 3d and space-like triangles in 4d). This 
summarizes the treatment and conventions of 112 111. The familiar form of a contribu- 



tion of a single bone b to the total curvature and therefore the action is 



A b S = vo\ume(b)5b 



where the volume of the bone is by definition real and positive and Sj, is the deficit 
angle around b. For the case of Euclidean angles this is conveniently illustrated by a 
two-dimensional example, where the bone is simply a point with volume 1. A pos- 
itive deficit angle 5 (Fig. [T5|a) implies a positive Gaussian curvature and therefore 
a positive contribution to the action, whereas an "excess" angle 5 contributes neg- 
atively (Fig. |l3|b). Lorentzian angles are also defined to be additive, in such a way 
that a complete rotation around a space-like bone gives 2n in the flat case. However, 
note that the angles are in general complex, and can become arbitrarily large in the 
vicinity of the light cone. In their contribution to the action, we have to distinguish 
between two cases. If the Lorentzian deficit angle 5 = rj is space-like (Fig. |i"4]a), it 
contributes as A b S = volume(6)?7b, just like in the Euclidean case. By contrast, if it 
is time-like (Fig. |T4lb), the deficit angle contributes with the opposite sign, that is, 
as A b S = —volume (b) r] b . Therefore, both a space-like defect and a time-like excess 
increase the action, where as a time-like defect or a space-like excess decrease it. 
The deficit angles in Sections 2.1 and 2.2 were calculated using 



cos 6 



(vi,v 2 ) 



sin 



a/ (Vi, Vi) (v 2 , v 2 ) - (vi, v 2 y 



9) 



(vi,vi)2{v 2 ,v 2 )* ' (vi,vi)i{v 2l v 2 )2 

for pairs of vectors V\,v 2 , and the flat Minkowskian scalar product (•, •). By defini 
tion, the square roots are positive imaginary for negative arguments. 
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Figure 14: Positive space-like (a) and time-like (b) Lorentzian deficit angles r\. 



Appendix 2: Link-reflection positivity in 2d 

In this appendix we will demonstrate the link-reflection positivity of the discrete 
model of two-dimensional Lorentzian quantum gravity introduced in [[7|. Recall 
that link reflection is the reflection at a plane of half-integer t. We choose it to lie 
at t = 1/2 and fix the boundary spatial lengths at the initial time t = —T + 1 and 
the final time t = T to 1-t+i = It — I- In order to describe a two-dimensional 
Lorentzian universe with these boundary conditions, we must not only specify the 
geometry of the spatial sections (simply given by l t , —T + 1 < t < T), but the 
connectivity of the entire 2d triangulation. 

A convenient way of parametrizing the connectivity that is symmetric with re- 
spect to incoming and outgoing triangles at any given slice of constant t is as follows. 
For any spatial slice at some integer time t, each of the l t spatial edges forms the base 
of one incoming and one outgoing triangle. The geometry of the adjoining sand- 
wiches is determined by how these triangles are glued together pairwise along their 
time-like edges. These gluing patterns correspond to distinct ordered partitions of 
the It triangles (either above or below t) into k groups, 1 < k < l t . We denote the 
partitions collectively by m(t) = {m r (t), r = 1, . . . , k} for the incoming triangles 
and by n{t) = {n r (t), r = 1, . . . , k 1 } for the outgoing triangles. The constraints on 

these variables are obviously Ylr=i m r(t) = Yl k r =i n r(t) — h an d a matching condi- 
tion for adjacent slices, namely, k'(t) — k(t + 1)Q In terms of these variables, the 

10 This parametrization is closely related to the description of 2d Lorentzian gravity in terms of 
"half-edges" (TsJ . 
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(unmarked) one-step propagator is given by 



fc>l n r ,m r >l, r=l,...,fc fe>l 



(90) 

where g = e x depends on the two-dimensional cosmological constant A. It is obvious 
from ( j90"D that the propagator depends symmetrically on the local variables m and 
n. The partition function for the entire system is 



II X^C'W. (91) 

t=-T+2 J t >l 

in accordance with (57). 

Following the discussion of link reflection in Sec. 6, we consider now functions 
F({m, n}) that depend only on the geometry above the reflecting plane at t — 1/2, 
i.e. on m(t), n(t), t > 1. The reflection (|52|) acts on the partition data according 
to 9i(m{tj) = n(—t + 1) and 9i(n(t)) = m(—t + 1). Without loss of generality, we 
can assume that F depends symmetrically on m(t) and n(t), Wt. Computing the 
expectation value (|54l) , we obtain 



((OiF)F) 

-l 

l ,h>l ' t=-T+2l t >l 

T-2 

(II E G 9 (k,lt + i)k + iG g (l T -i,lT)F) 
t=l l t+1 >l 

= & e «. e 1 (' ° : i) G; : i) 

/c>i 

Since the last expression is a sum of positive terms, we have hereby shown that 2d 
Lorentzian gravity is link-reflection positive. 
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